Hybrid theory of electron-hydrogenic systems elastic scattering 


A. K. Bhatia 

NASA/Goddard Space Flight Center, Greenbelt, Maryland 20771 USA 


Introduction 

Accurate electron-hydrogen and electron-hydrogenic cross sections are 
required to interpret fusion experiments, laboratory plasma physics and 
properties of the solar and astrophysical plasmas. 

The advantage of working with hydrogenic targets is that the wave function 
is known exactly which allows us to test any new theories on scattering of 
the simplest three-body systems. The incident electron produces a distortion 
of the target orbital and various methods have been used to take into account 
this distortion producing a long-range potential. Among them is the method 
of polarized orbitals [1]. In spite of the fact that this method includes the 
essential physics by modifying the target wave function suitably in the 
presence of the incident electron, the results obtained for the phase shifts do 
not obey any bound principles. In spite of this shortcoming, this method has 
been used extensively on electron scattering from atoms [1] and molecules 
[ 2 ]. 

In spite of singularities in his calculations, Schwartz [3] used the Kohn 
variational principle to obtain accurate phase shifts for electron-hydrogen 
scattering. More recently, scattering from hydrogenic systems has been 
carried out by using the Feshbach projection operator formalism [4]. In these 
calculations [5, 6, 7], the short-range correlations could be included 
explicitly by introducing separate correlation functions and then 
amalgamating them into the scattering equation via an optical potential, thus 
replacing the many- particle Schrodinger equation with a single-particle 
equation. The results obtained for the phase shifts are very accurate and they 
have an important property namely that they have the rigorous lower bounds 
to the exact phase shifts. That is the exact phase shifts are always higher than 
the calculated phase shifts. In this way, we obtained accurate -S-wave and P- 
wave phase shifts for electron-hydrogenic systems. This approach also helps 
us to get very accurate scattering wave functions for S and P states in He + 
ion and H atoms. Very accurate cross sections for photoionization of He and 
photodetachment cross sections of H' ion were obtained along with radiative 



attachment cross sections [8], A detailed account of these methods and 
calculations has been given in [9], 

In these calculations [5, 6, 7, 8], the short-range correlations could be 
included explicitly but not the long-range correlations at the same time. In 
this article, we present an approach which includes the long-range 
correlations as well as short-range correlations explicitly at the same time 
[10]. Also, we do not use the projection operator formalism of Feshbach [4], 
Let us write the wave function for an angular momentum L in the form 

+ ( 1 ) 

r \ 


The summation index X takes values from 1 to the number of terms N in the 
wave function O, the correlation function. The target function is given by 



Z is the nuclear charge of the target. The (±) above in Eq. (1) refers to 
singlet (upper sign) or triplet (lower sign) scattering, respectively. The first 
term represents the static approximation. Beyond the first term are the terms 
giving rise to the exchange approximation [11], The function is the 
correlation function which for arbitrary L can be written in terms of 
symmetric Euler angles [12]: 


®z=Z[/f’ +1 (n 


* )D1 M (A <p, w) + ft ' E, , r 2 , r n )T>r 1 (</>, <p, W )] 


(3) 


The D functions have been called rotational harmonics and 9, cp and \|/ are 
the Euler angles. The/’s above are the generalized ‘radial functions’ which 
depend on the three residual coordinates which are required to define the 
two vectors iq and r 2 . The radial functions have to be defined for each 
angular momentum L [12]. 

The Euler angles under exchange are given by 


& ^ n - 6, (p ^ n + (p,^ ^ 2n 


(4) 



The Hamiltonian in Rydberg units can be written as 


H = -V 2 -V 2 2 


2 Z 

r i 



( 5 ) 


and the total energy is given by E=k 2 - Z 2 , where k 2 is the kinetic energy of 
the incident electron. 


Taking N=\, for illustration, calculate the functional (arising in the Kohn 
variational principle) 


I = <\|/ L (r i , r 2 )|H-E| \|/ L (r i,r 2 )>, (6) 

which can be written as 

I = A + C X B + C 2 D, (7) 

where 

A=<[cp 0 (r 2 )u(ri)±(l <^2)]|H-E|cp 0 (r 2 )u(ri)±(l <->2)]) (8) 

and 

B=2< ® f (t\,r 2 ) \ H - E \ (P 0 (r 2 )u(r, ) ± (1 2 )> (9) 

= 4<V l (r l )u(r l )> 

In the above equation, 

E(h) =< \H-E\ <p 0 ( r 2 ) > . (10) 

In Eq. (7) 

D=<®?(r l ,r 2 )\H-E\<t>V(r l ,r 2 )>=s l -E, (11) 

where Si is the expectation value of H 

\H\0>f >=£!. (12) 



The wave function is normalized to unity. Now we determine the 
eigenvector Ci. The variation with respect to Ci in Eq. (7) is 


dl 

dC x 


= 0, 


Implies 


(13) 


B + 2C X D = 0. 


(14) 


This gives 

c B _ 2<&"(r l ,r 1 )\H -E\(p {2 {r 1 )u L {T\)> 

1 2D E-e l 


2 < V\( r i)u L (i\) > 

E -Sj 


(15) 


Now the correlation term in Eq.(l) is known. Thus the variational principle 
reduces to 


<cp 0 (r 2 ) \H-E\ \)/l> = 0 , (16) 

which can be simplified to 

< <P 0 (r 2 ) I H - E I [<p 0 (r 2 )u(T\ ) ± (1 2)] > +C,f (F t ) = 0 . (17) 

Where Vi(ri) is given in Eq. (10). Substitution of Ci from Eq. (15) into Eq. 
(17) gives the equation 


< <P 0 (r 2 ) | H -E | [(p 0 (f 2 )u(r t ) ± (1 2)] > + 2V ^^ < V ^ ^ 1 > = 0 

E - £ j 


(18) 


Now the resulting equation for u L { r), letting ri=r, can be written in the form 


+ 1) + ^ ±Va _ v k 2 K =0 , 
dr" r" 


(19) 



where 


tf E-e i 


( 20 ) 


The above equation also holds when X ranges from 1 to N, the number of 
terms in the correlation function (3). The potentials V d is the direct potential 
and V ex is the exchange potential. The optical potential has been obtained 
without the use of the Feshbach projection operator formalism and hence is 
independent of the projection operators P and Q [4]. The phase shifts 
obtained by two approaches are the same within the numerical accuracy. 


Optical potential with polarization 


To take into account the long range correlations, we replace ( r 2 ) in Eq.(l) 
by 


<b pol (r,,r 2 ) = <f>(r 2 ) 


s{i\,r 2 ) u is ^ p (r 2 ) cos (0 l2 ) 
r, 2 r 2 


( 21 ) 


where 


is—>p 


(r 2 ) = e 


-Zr 2 


Zr 

^__ 2 _ 


+ K ), 


( 22 ) 


which is the dipole part of the polarized orbital and 0i 2 is the angle between 
r! and r 2 . Temkin [1] introduces the function 


e(r 1 ,r 2 )=l, r!>r 2 

= 0, r, < r 2 , (23) 

which ensures that the polarization takes place when the scattered electron is 
outside the orbital electron. Now we want the polarization function in Eq. 
(21) to be valid through out the range, rather than only just for r! > r 2 . We, 
therefore, replace the function s(r b r 2 ) by a cut off function and in order to 
avoid discontinuity at r 2 we choose it of the form 



Xp {r x ) = { l-e^y, 


( 24 ) 


where the exponent n > 3. The nonlinear parameter ft, which is function of 
k, can be used to optimize the results. This function vanishes at large values 
of r 1? it also contributes to the short-range correlations in addition to those 
obtained from the correlation function <f> L . Another form of the cutoff has 
been proposed by Shertzer and Temkin [13] 


Xst = d 3 r 2 0 o (r 2 ) ^ 

N J r, r, 


(25) 


where N denotes the normalization and is given by 

, T f ,3 \ U \s->p ( r 2 ) 1 

N = \d r 2 <f> () (r 2 ) 

J r r 


(26) 


In the above equations the integrations are from 0 to infinity. This gives 

fe (r 1 ) = l-e- 2Z 'T^(Zr 1 ) 4 +|(Zr 1 ) 3 +2(Zr 1 ) 2 +2Zr 1+ l], (27) 

which guarantees that ^r^/rf approaches 0 when ri approaches 0. Now the 
scattering wave function can be written as 


T L (r 1 ,r 2 ) = 




Y lo (n, )® r 2 )± (1 « 2) + C, © J (r, ,r, ), 


(28) 


where the target function with polarization is given by 


(b pol (r ] ,r 2 ) = </>(r 2 ) 


z(n) u is—>p (r 2 ) cos(0 n ) 
K r 2 


(29) 


The cutoff function % can be of either form given in Eq. (24) or Eq. (27). 

We arrive at the same form of Eq. (18) by replacing cp 0 by ® po! (r x , r 2 ) in Eqs. 
(6)-(20). We restrict ourselves to L = 0 and we can write the final equation in 
the form 



(30) 


[D(r)^j + k 2 +V, +V-' ±(T„ + V*)-V*Mr) = 0 


The various quantities are given below: 


D(r) = \ + ^-(^X) 2 . 
K 8 Z 6 r 2 


The direct potential is given by 


V d = 2(Z ^ + 2e~ 1Zr (Z + -) 


1 

r 


and 


v f° l = ( x i +x,) + x 2 -^~. 

dr 






43 ,x(r\ , 2a(r) . x 4.5^(r) 2 8 ,,, r x(r) 


8Z r (Zr) 


d(r)(^-y 

(Zrf 3 Z r 2 


5, (r) = 4-[-2Z^(r) + 2Z - - + e- 2Zr (- |z V - 1 Z V 2 + 2Z + -)], 


{ \ T7 0 / ^ 4Z 6 -22t/-4 v 5 , tv 4 4Z 8Z 6 

B 2 (r) = -2ZB 1 (r) r + — + e (-Z r + 2Z 


The polarizability of the target as a function of r is given by 


a(r)/Z 4 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


where 



(40) 


a(r) = 4.5-|e“ 2Zr (|(Zr) 4 +7.5(Zr) 3 + 13.5(Zr) 2 + 13.5(Zr) + 6.75). 


and 


d(r) = 


129 1 18 1 2 3 4 27 

— -e 2zr (—r 4 + r 

32 Z 5 r Z 7 r 3 16 16Z 



135 975 1281 1 36 1 18 1 

H r f “I T "I t 1 7 — H t — r) 

8Z 3 32Z 4 32Z 5 r Z 6 r 2 Z 7 r 3 


(41) 


*3 -> r infinity 

In the above expression 9/(2Z 4 ) is the dipole polarizability of the target with 
nuclear charge Z. The direct polarization potential has also terms 
proportional to 1/r 2 and 1/r 3 which go to zero for r ->0. 

For completeness, we give other terms in the Eq. (30). The exchange terms 
are given by 

V ex u(r) = 4 Z 3 e~ Zl [(/c 2 + Z 2 )rj dxe~ Zx xu(x)~ 2 dxe~ Zx u(x)-2rj dxe Zx u(x)j, (42) 

Vf x ol ii(r) = J 2 u(r) + J 3 u(r) + J 4 u(r), (43) 

where 


J 2 u(r)= |S >j,( \ -f dxe Zx %(x)u(x) + r 2 f dxe z 4(i)^], 
3 r r Jo Jr x 


(44) 


J 3 u(r) = 


8Ze z ' Z (r) 1 

3 r 2 r 


[~ J 0 dxuu-^p (x)xu(x) + r 2 £ dxu u ^ p (x)xu(x) + r 2 J dxu ls ^ p (x) -XI ] 


u(x) 1 

x 2 

(45) 


and 


8 Ze- Zr Z (r) t 


[J 0 dxUu-> P ( x )^T-Z(x)\ + G 2 u(r) + G 3 u(r). 

X 


J 4 u(r ) = 


3 


(46) 



Now we will write below the expressions for G 2 u and G 3 u. 


G 2 u(r) = 


4 Z(r)u u ^ p (r) 


3Zr 


1 f 00 V 

- 1 tlxe z ' [ 

Jo 


3 ,.3 


Zx 


-2 Z 2 x 2 -Zx + 2 + (— + x z )(— ~^ + k z )y 


. Zx 3 , 2w 2Z 2 _ ,, 2vi *(*)«(*) 




(47) 


%Z( r ) u u^> n ( r ) 1 r u \ r® 

= — ^ [- f Z(x)u(x) + f dx — z( x ) u (x) 

8Zr r Jo x Jr x 


i + ~Gu(r)\ 

(48) 


Gu(r) = -^ £ dxu u ^ p {x)%{x)u{x) + r 2 £ dx U ^ P S %{x)u{x) 


(48a) 


The polarization term in Eq. (21) gives rise to terms in equations (33 - 41) 
and equations (43-48). These terms make the calculations much more 
intricate compared to the calculations without these term. 

The optical appearing in Eq. (30) is given by 


Nco 


Kfp'uin ) = Yu 


<Y m (n i )^ ol (? l ,r 2 )\H-E 




x ® ( 0 S) 


H-E I > 


E-s „ 


(49) 


It should be noted that in the above expression, we have Hy instead of Hf, 
given in Eq. (1). The function H',', is the wave function without the 
correlation terms in Eq. (1). This is very clear in the projection operator - 
formalism of Feshbach [4] used in [5, 6, 7, 8] to obtain the scattering 
equations. There the function Hy is written as P V F 0 + Q V F 0 and the function 
P y Y () only occurs on the right hand side of the formulation. It should be 
pointed out that P+Q= 1 . Now, we have given most of the quantities 
occurring in Eq. (30). 

For -S-waves (i.e., L— 0), /^constant in Eq. (3) and the correlation function 
is a function of “radial” coordinates. We take O 0 of Hylleraas form 



(50) 


N m 


(\i l=n (? l j 2 ) = Y J c , mi i^ 


yt i 


-S Tj I 


12 


- (1 <— 2 )], 


Imn 


The sum includes all values of l , m and n such that l+m+n=(o and co is equal 
to 0, 1, 2,3 . .. 8 . The total number of terms N m depends on spin and whether 

y=5 or not. There are 95 terms for l S state and 84 terms for 3 S state for N (0 = 8 
in these calculations. Solutions of the integrodifferential equations have been 
carried by the noniterative method [see Appendix A] and by using 
quadruple precision. We present results for the S and S states in Table I 
and II in three approximations: 

A. Phase shifts calculated using optical potential obtained by using the 
projection operator formalism [4], 

B. Phase shifts calculated using optical potential obtained without using the 
projection operator formalism [4], Eq. (19). 

C. Phase shifts calculated using optical potential without using the projection 
operator formalism but including the effects of polarization, Eq. (30). 


Table I. Phase shifts (radians) of l S state for various k. 



Ref. 5 


C 

C 

Ref. 3 

Ref. 15 

k 

A 

B 

XST 

Xb 

h Schwartz 

PSSB 

0.0 a 


6.05327 

6.00092 

5.99567 

5.965 


0.1 

2.55358 

2.55158 

2.55372 

2.55370 

2.553 

2.550 

0.2 

2.06678 

2.06644 

2.06699 

2.06717 

2.0673 

2.062 

0.3 

1.69816 

1.69640 

1.69853 

1.69684 

1.6964 

1.691 

0.4 

1.41540 

1.41783 

1.41561 

1.41554 

1.4146 

1.410 

0.5 

1.20094 

1.20084 

1.20112 

1.20195 

1.202 

1.196 

0.6 

1.04083 

1.04074 

1.04110 

1.04191 

1.041 

1.035 

0.7 

0.93111 

0.93105 

0.93094 

0.93084 

0.930 

0.925 

0.8 

0.88718 

0.88717 

0.88768 

0.88802 

0.886 



a k=0 results represent scattering 


engths. 


The scattering length a is defined by 


lim^o k cot rj = -\ / a 


(51) 



It should be noted that the phase shifts given above have the rigorous lower 
bounds to the exact phase shifts while the scattering lengths have the 
rigorous upper bounds to the exact scattering lengths. 

Temkin [14] has shown that the scattering length is significantly affected by 
the long-range polarization. The scattering length obtained at a short 
distance R can be corrected by the following formula 

a = a(R)-a(j- — + -—), (52) 

where a is the true scattering length. The scattering length a(R) for l S state 
obtained using xp is 5.99567 at 7?=146.0782 for N=95 terms. Using the 
formula (52), we get «=5.9661 1 which agrees well with the value 
a =5.965 ± 0.0003 given by Schwartz [3], He also included long-range 
polarization in his zero-energy calculation. The scattering lengths are 
rigorous upper bounds to the exact scattering results. 

We find that the phase shifts obtained with the Feshbach formalism [4] are 
fairly close to the present results. That indicates that the short-range 
correlations do try to assimilate the effect of the long-range correlations to 
some extent. The present phase shifts have rigorous lower bounds to the 
exact results. We compare our results with those of Schwartz [3], His 
calculation was carried out using the Kohn variational principle. But his 
calculations do not provide any bounds on the phase shifts though the Kohn 
variational principle does provide a bound on the scattering length. 

In Table II, we give the results for the S state. The S scattering length for 
W=84 with x P , [3=0.6 and n= 4 is a(R)= 1.7681542 at/?=349.0831. The 
corrected value, using formula (52), is 1.76815 is a little lower compared to 
the Schwartz value 1.7686 ±0.0002. The formula (52) certainly gives a good 
agreement with the Schwartz value. Lower scattering lengths are consistent 
with the rigorous upper bounds, as state earlier. 

In tables I and II, we have given rj SSB obtained by Scholz, Scott and Burke 
[15] using the well known 7?-matrix formulation. In these calculations, all 
the correlations are included in the inner region of a certain radius R. The 



functions in the two regions are then matched at this radius R. Their results 
are fairly close to the presently calculated phase shifts. 


Table II. Phase shifts (radians) of S state for various k. 



Ref. 5 


C 

C 

Ref. 3 

Ref. 15 

k 

A 

B 

XST 

X* 

h Schwartz 

PSSB 

0.0 a 


1.81644 

1.78467 

1.78154 

1.7686 


0.1 

2.93853 

2.93850 

2.93856 

2.93856 

2.9388 

2.939 

0.2 

2.71741 

2.71740 

2.71751 

2.71751 

2.7171 

2.717 

0.3 

2.49975 

2.49956 

2.49987 

2.49987 

2.4996 

2.500 

0.4 

2.29408 

2.29394 

2.29465 

2.29457 

2.2938 

2.294 

0.5 

2.10454 

2.10414 

2.10544 

2.10574 

2.1046 

2.105 

0.6 

1.93272 

1.93280 

1.93322 

1.93336 

1.9329 

1.933 

0.7 

1.77950 

1.78010 

1.77998 

1.77998 

1.7797 

1.780 

0.8 

1.64379 

1.64408 

1.64425 

1.64444 

1.643 



a k=0 results represent scattering 


engths. 


We see that the hybrid theory developed here includes both the short-range 
and long-range forces. The present development is significant even when the 
present results do not differ significantly from those obtained using the 
Feshbach formalism. 


Electron-He + Scattering 

The above approach has been applied to the scattering of electrons from 
helium ions. Now the phase shifts are given by 

lim,^ u(r ) oc sin[£r + — — -ln(2£r) + argT(l - — — ) + 77]. 

k k 

(53) 

Phase shifts obtained for e-He in S and S states are given in Tables III and 
IV. 



They are given at low energies as well and have been compared with the 
results obtained using the Feshbach projection operator formalism [6]. They 
are also compared with those obtained in the close-coupling approximation 
[16] and also obtained using the Harris-Kohn method [17]. 


Table III. Comparison of the presently obtained *S phase shifts in e-He + 
scattering with other calculations: (OP) Ref. 6, (CC) Ref. 16 and Harris- 
Kohn Ref. 17. 


k 

Hybrid 

(OP) 

(CC) 

Harris-Kohn 

0.1 

0.43808 




0.2 

0.43550 




0.3 

0.43142 



0.4300 

0.4 

0.42608 

0.42601 


0.4228 

0.5 

0.41974 

0.41964 


0.4078 

0.6 

0.41265 

0.41278 

0.4111 

0.4086 

0.7 

0.40568 

0.40561 

0.4046 

0.4024 

0.8 

0.39865 

0.39857 

0.3974 

0.3968 

0.9 

0.39213 

0.39202 

0.3906 

0.3893 

1.0 

0.38644 

0.38634 

0.3850 

0.3836 

1.1 

0.38200 

0.38187 

0.3805 

0.3794 

1.2 

0.37914 

0.37899 

0.3780 

0.3741 

1.3 

0.37846 

0.37832 

0.3744 

0.3721 

1.4 

0.38158 

0.38560 


0.3786 

1.5 

0.39802 



0.4014 

1.6 

0.34480 





1 3 

Phase shifts are higher for triplet states and both S and S phase shifts 
decrease with increasing incident energy. It is seen that phase shifts start 
rising after k=\A, indicating the existence of a resonance which we will 
discuss below. 




Table IV. Comparison of the presently obtained J S phase shifts in e-He + 
scattering with other calculations: (OP) Ref. 6, (CC) Ref. 16 and Harris- 
Kohn Ref. 17. 


k 

Hybrid 

(OP) 

(CC) 

Harris-Kohn 

0.1 

0.93065 




0.2 

0.92704 



0.9270 

0.3 

0.92114 



0.9210 

0.4 

0.91302 

0.91300 


0.9128 

0.5 

0.90282 

0.90275 

0.9019 

0.9025 

0.6 

0.89057 

0.89050 

0.8910 

0.8902 

0.7 

0.87645 

0.87640 

0.8777 

0.8762 

0.8 

0.86066 

0.86069 

0.8617 

0.8605 

0.9 

0.84366 

0.84356 

0.8440 

0.8435 

1.0 

0.82536 

0.82531 

0.8253 

0.8251 

1.1 

0.80636 

0.80625 

0.8062 

0.8062 

1.2 

0.78677 

0.78666 

0.7868 

0.7865 

1.3 

0.76696 

0.76684 

0.7672 

0.7665 

1.4 

0.74708 

0.74697 


0.7466 

1.5 

0.72746 



0.7274 

1.6 

0.70815 



0.7095 


Electron-Li ++ scattering 

A similar calculation has been carried out for e-Li ++ scattering. As in the 
case of He + , there is a long-range Coulomb potential in addition to the long- 
range potential due to polarization of the target. 

Gien [18] has carried out calculations by using Harris-Nesbet method and 
has used various combinations of target states in his calculations. We have 
given results which are labeled (E4S) in his paper. He used Is, 2s, 2 p states 
of the target along with a 2 p pseudostate [19] and correlations. The 
pseudostate is given by 

2^ = Z V * 7I [0.340rV 0 ' 5Z ' - 0.966r V Zr (l + 0.5Zr)]. 


(54) 




This state has been formed as a linear combination of target 2p state and 
U| S _> P with a constant factor such that it is normalized and orthogonal to the 
2/7 state. The inclusion of the pseudostate gives the exact polarizability of the 
target. 


1 3 

Table V. Comparison of the presently calculated S and S phase shifts in e- 


Li ++ scattering with those obtained by Gien [18]. 


k 

^(Hybrid) 

! S(Gien) 

3 S (Hybrid) 

3 S(Gien) 

0.1 

0.23188 


0.56084 


0.2 

0.23176 


0.56020 


0.3 

0.23148 


0.55869 


0.4 

0.23109 


0.55678 


0.5 

0.23064 

0.2273 

0.55435 

0.5526 

0.6 

0.23012 

0.2264 

0.55142 

0.5499 

0.7 

0.22960 

0.2265 

0.54799 

0.5467 

0.8 

0.22906 

0.2272 

0.54413 

0.5430 

0.9 

0.22855 

0.2277 

0.53925 

0.5390 

1.0 

0.22807 

0.2275 

0.53514 

0.5345 

1.1 

0.22769 

0.2262 

0.53000 

0.5296 

1.2 

0.22740 

0.2250 

0.52456 

0.5244 

1.3 

0.22724 

0.2258 

0.51880 

0.5189 

1.4 

0.22724 

0.2328 

0.51276 

0.5131 

1.5 

0.22742 

0.2521 

0.50646 

0.5069 

1.6 

0.22782 


0.49997 

0.5005 


Gien [17, 18] obtained phase shifts at irregular energy points, obtained from 
the diagonalization of the Hamiltonian His results could be fitted to 

>l = A + BE + CE 2 +DE i +Fexp[aE], (55) 

In the above equation E=k 2 is the incident energy and A, B, C, Z), E, E and a 
are the fitting parameters. 



Resonances 


The first time, a resonance was calculated, was in 1962 by Burke and Schey 
[20], This was in electron-hydrogen scattering. They used the close-coupling 
approximation and found that the resonance was centered at 9.61 eV. 

Resonances have been known by many names: autoionization states, doubly 
excited states and Feshbach resonances. It is observed that there is a rapid 
change of phase shift in the resonance region. This phase shift can be fitted 
to the Breit-Wigner form to obtain the resonance parameters, 


n C aicX E ) = n 0 +h£ + tan 


0.5T 

(Er-E)’ 


(56) 


2 

where E=k is the incident energy, rj calc are the calculated phase shifts, r| 0 , 

A, T, and E R are the fitting parameters. E R and T are the resonance position 
and resonance width. We would like to calculate the resonance position and 
width of the lowest 'S resonance in He (below the n — 2 threshold of He). In 
Table VI, we give phase shifts for 70 terms in the correlation function in the 
resonance region. 


Table VI. Phase shifts in the resonance region for the lowest 'S resonance in 
He and Li + . 


Target 

k 

h 

Target 

k 

h 

He + 

1.555 

0.5483195 

l r 

2.282 

-0.090487 


1.558 

0.6398967 


2.2825 

-0.058190 


1.560 

0.7803777 


2.283 

-0.031793 


1.562 

1.177147 


2.2835 

-0.009842 


1.5634 

1.938293 


2.284 

0.008626 


1.56345 

1.972649 


2.2845 

0.0.02451 


1.565 

2.776815 


2.285 

0.038199 


1.5655 

-3.372919 





1.566 

-3.276240 





1.567 

-3.148991 





We calculate the square of the difference between the left and right side of 
Eq. (57). The difference is being calculated at the energies given in Table 




VI. Minimize the sum of the squares for various fitting parameters. We find 
r| 0 =0.3761, A=0, 2t R =2.4426 Ry with respect to He + and r=0. 00906 Ry when 
the minimized sum is 1.57xl0' 7 . This gives resonance position 2i R =57.8481 
eV with respect to the ground state of He and width T=0. 1233 eV. This 
agrees very well with 2^=57.8435 eV and T=0.125 e V, obtained using the 
Feshbach formalism [21], These resonance parameters, calculated now, 
include the contribution from the long-range and short-range correlations. 
But these parameters do not have any rigorous bonds although the calculated 
phase shifts do have rigorous lower bonds. This is the disadvantage of fitting 
to the Breit-Wigner form Eq. (56). 

A similar calculation has been carried out for the lowest 'S state in Li + . The 
phase shits in the resonance region for 70 terms in the correlation function 
are given in Table VI. By fitting, we obtain 2i R =70.5904 eV with respect to 
ground state of Li ++ and T=0. 1657 eV which compares very well with 
2i R =70.5837 eV and T=0.157 eV [22]. In the Feshbach formalism, the 
resonance parameters, position and width, are defined as 

E r = <0|QHQ|T» + A = 8 R +A (57) 

r=2£<\|/|PHQ|T» (58) 


It should be noted in the projection operator formalism of Feshbach, the 
contribution A to the position of the resonance from the continuum and 
nearby resonance has to be calculated separately. This contribution is 
already included in the present calculation which is a great advantage. 


Conclusions 

We have developed a method in which the short-range and long-range 
correlations can be included at the same time in the scattering equations. The 
phase shifts have rigorous lower bounds and the scattering lengths have 
rigorous upper bounds. The phase shifts in the resonance region can be used 
to calculate very accurately the resonance parameters. 



Appendix A 


We will briefly describe the non-iterative method for solving 
integrodifferential equations [23], Consider the equation for the scattering 
function u(r) given by 

[^j + V(r) + k 2 ]u(r) = g(r)\ Q f(x)u(x)dx (Al) 

Let 

u(r)=u 0 (r) + Cui(r) (A2) 

The constant C represents the definite integral in the Eq. (Al). 

Substituting u(r) in Eq. (Al), we get two equations 


d_ 

dr 2 


+ V (r) + k 2 ]u 0 (r) = 0 


(A3) 


d_ 
dr 2 


+ V(r) + k 2 ]u 1 (r) = g(r) 


(A4) 


These two equations can be solved easily. The substitution of (A2) in C 
gives 

C = [ f(x)u 0 (x)dx+c\ f(x)u l (x)dx = I 0 +CI l (A5) 

JO JO 

Having calculated u 0 and Ui, I 0 and f can be calculated. From the above 
equation, we can now solve for C. We get 


C = 


(i-A) 


(A6) 


The Eq. (Al) can now be written as 

+ V(r) + k 2 ]u(r) = Cg(r ) 
dr" 


(Al) 



Now the right-hand side is known and this equation can be solved for u(r). 
This method can be generalized to any number of definite integrals. 
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